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Abstract 

Based on the Stochastic Inclined Rods Model (SIRM) proposed by 
H. Matsuura and M. Nakano, we study the microscopic motion of actin 
myosin system including the motion of the G-actin. Our model is com- 
posed of an inclined spring (rod), a myosin head, a myosin filament and 
G-actins. We discuss the stochastic resonance between the myosin head 
and random noise. The results of calculation show that the model can 
convert the random motion to one directional motion, and the myosin 
head works as a resonator of random noise which absorbs the energy 
through the stochastic resonance. The intermolecular potential between 
the myosin head and G-actin and the inclined rod play a key role for the 
muscle's motion. The energy consumed by the motor is directly supplied 
from the surroundings (i.e., the thermal motions of water molecules). 
Keywords : SIRM, Intermolecular potential, Actin myosin system. 



1 Introduction 

Biological systems are among the most challenging subjects for theoretical physi- 
cists as well as experimentalists or simulationists. More and more biologists 
and theoretical biophysicists have focused on studying the mechanism of living 
systems and mincing the biological structure to invent subtle artificial instru- 
ments. Nowadays it is exciting to study a problem how motor protein gets a 
energy of movement and what makes myosin slide in one direction along the 
filaments. Based on different system, three different families of motor proteins 
have been identified PP |Ej '■ Kinesis and dyneins move along tubulin filaments 
0], myosin moves along actin filaments The motion mechanism of these 

motors can be described in general as follows pQ [7] [8] : firstly the motor protein 
binds Adenosinetriphoshpate (ATP), then hydrolyzes the bound ATP and ab- 
sorbs the energy, subsequently, it releases the products Adenosinediphosphatc 
(ADP) and Phosphate (P). So the motors move constantly within chemical cy- 
cle. The molecular motors play a key role in transducing chemical energy into 
mechanical work at a molecular scale. 

In 1954 H. F. Huxley, H. E. Huxley and R. M. Simmons HHdOHm proposed 
the rotation cross bridge model to study muscular mobility. In these earliest 
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theoretical descriptions, the asymmetry of the system was introduced via asym- 
metric transition rates and specific conformational changes were also discussed. 
This model has been widely spread in past two decades and has become a 
standard model even in college textbooks, but extensive experimentations have 
failed to find evidence for a large conformational change in myosin head during 
the force generation process |12| . So it is worth while to reconsider an old idea, 
proposed by A. F. Huxley in 1957, the motor protein uses thermal energy that 
exists at physiological temperature to perform work |13| . such as the thermal 
ratchet model ^1] ^Hl, the isothermal ratchet model [16| . 

Recently, there have been many reports that discuss the relation between 
the muscle's movement and Stochastic Resonance (SR) or thermal noise JZ] Q2| 
J^IEni- These explanations, however, somewhat abstract from the viewpoint 
of biology, and these models are lacking in dynamic descriptions of actin myosin 
interaction. 

The newly resonance models for muscle's motion was named the Stochas- 
tic Inclined Rods Model (SIRM) proposed by H. Matsuura and M. Nakano 
[2"T| |22 [221 EU • The energy of the motion was supplied from the random noise 
and the system always moves to one direction by using of stochastic resonance. 
The movement of the system does not break the second law of the thermody- 
namics, because the actin myosin system is open to the surroundings and the 
energy flows in from their surroundings. The SIRM presents a perfectly sliding 
mechanism for the actin myosin system. 

The main aim of our paper is to study the motion of actin myosin system in- 
cluding the G-actin's motion. We discuss the stochastic resonance of the myosin 
head and give the numerical results for the motion of the system. We find that 
there is a relative sliding motion between the G-actin and the myosin owing to 
the intermolecular potential and the structure of the inclined rod. 



2 Model and Formalism 

The structure of the actin myosin system is shown in Fig. l[2I] together with 
structure of the muscle. It consists of the chain structure of the actin of length 
around 100 nm and a bundle of hundreds of myosin which has an inclined rod 
and a head of 0.45 million Dalton. The whole system exists in a water solution 
and the interaction with water molecules cannot be neglected. 



Fig. 1 

In order to explain the contraction of the muscle, namely solve the motion 
of actin myosin system, we simulate it by a mechanical model as shown in Fig. 
2. We call this model as stochastic inclined rod model (SIRM). The motion 
mechanism of the SIRM is as follows: firstly, ATP is hydrolyzed and releases 
the energy which makes the random noise interacts with the myosin head and 
myosin head obtains the energy from random noise then the head vibrates, it 
collides with a G-actin and obliquely kicks the G-actin sphere, because the di- 
rection of vibration is inclined against the line of the actin fibers and myosin 
molecules obtain the propellant force along the direction of the fibers. In this 
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way, the filament can move to one direction while the G-actin moves to oppo- 
site direction. The actin myosin system can move to one direction owing to the 
structure of inclined rods of myosin and the intermolecular potential between 
the myosin head and the G-acin. 



Fig. 2 

According to Fig. 2 we can construct the equations of the motion for the head 
of myosin. 

With regard to the center of gravity (x2, yi) of the filament, we set a similar 
equation: 

The variable y 2 is fixed since the myosin filament does not significantly move 
along the y-direction compared to the re-direction. 
The equation of the G-actin is shown as follow: 



9 dt 2 ~ dx g 7 m ■ 



Where F x (t), F y (t),F X2 (t), F x (t) are fluctuation of the thermal noise, m, M, M g 
are the masses of myosin head, myosin filament and the G-actin, respectively, 
a, /3, rj, 7 are viscous constants, and the fluctuation forces follow fluctuation- 
dissipative relation |2"5] |2"f)] [?7j : 

< F a {t) >= 0. (5) 

< F a (t)F b (s) >= 2k B T(6 a , b 6(t - s). (6) 

where a, b = x, y, Xi, x g ; Q = a, (3, rj, 7; ks is Boltzmann constant. T is absolute 
temperature, t, s are time. 

As for the intermolecular potential between the myosin head and G-acitn we 
adopt Jcnnard-Jone potential: 



U a = J2(P^ 12 -Q^ 6 )U. (7) 



Where n = \/ (x - Xi) 2 + (y - y 4 ) 2 - R; x { = x g + iR; y t = y g ; (x g ,y g ) 
stands the coordinates of the center of the G-actin and R is its radius, p, q, U 
are the parameters of the potential. 

The potential of the myosin rod is approximately expressed as follows: 



U s = K t exp{-^{x-x 2 )i + (y-y 2 ) 2 + L) 
+Ki exp(^(x-x 2 ) 2 + (y-y 2 ) 2 - L) 

+ l -K e {6-e»f. (8) 
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Where Kg,Ki are parameters of potential of the spring and L is nature 
length of the spring, tan(0) = X2 ~ x , 6, 6>o are current angle and initial angle be- 
tween the rod and level axis. (2:2,2/2) * s ^ ne cen ^ eT °f gravity of myosin bundle. 



3 stochastic resonance 

We discuss the stochastic resonance between the myosin head and the thermal 
noise to show that SIRM moves to one direction with accepting energy from 
thermal noise. Because the above equations are so complex and highly non- 
linear that we cannot obtain a strict analytical solution, we only discuss the 
motion of the myosin head and set up the following simple equation for the 
head, where the j-th G-actin is the nearest one to the head: 

x dx 
m— = 12pU(x - Xj)- 13 - 6qU(x - XjY 1 - 2(x -x 2 + L)K + F x (t) - a—. 

(9) 

Because the fluctuation F(t) can be given in the following form |22| : 



F a {t) = ^M^^expW). (10) 

3 

< F a (t)F b (s) >= ^-^Y^S^expiirts-QH). (11) 

3,1 

< (CPb >= S a,b, < Pi* > = < p\ >= 0. (12) 

If we assume that the thermal noise contains all of frequency, one of the 
noise frequency is expressed as follows: 

EM _> Zexp(ifit). (13) 

TO TO 

Thus Eq.(9) can be given: 

-± + 2^ + u 2 i = J- exp(tfii) - eai 2 - sfc 3 - (14) 
ot z at to 

Where £ = x - x 2 + L, = Xj - x 2 + L, 2/i = 

2 _ 2K 156p[/ 42qU _ 21M P U 336qU _ 32760pU 3024qU 
W ° ~ ~ + ^j r ~^f ,£a ~ 2mtf _ "2^eT' £/ ^ _ 3!to£ 16 ~ 3!to£)° ' 

' / ^ . (15) 

Here we use Taylor expansion for the intermolecular force of the G-actin and 

ignore the higher order terms. 

When e = and /j 2 > ujq we can get the solution of the Eq. (14) as follows: 



£ = aoe - "' cos(ut + ti ) + ae l{nt - s) . (16) 
- P 2 : a = - 

Ti 

determined by initial conditions. 



Where w = ^J^-p\ a = m ^ {ujl J 2)2+i ^ tan 6 = a 0l i9 are 
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If e ^ 0, we can use a effective frequency for the Eq. (14) 

3ea 2 

So we can rewrite the Eq.(14) as follow: 



= s-r. (17) 



d J + 2,§ + ^ = l cM m- (18) 

So the solution of the Eq. (18) should have the same form as that of Eq.(14), 
but the amplitude a should be satisfied the following relation: 

a 2 = l™i (19) 

If we assume the frequency Q = luq + A, we can obtain: 

K 2 -tt 2 ) 2 =4^ (!^-A) 2 . (20) 
We can rewrite the Eq. (19) as follow: 

a[(a-A) 2 + (^)] = F. (21) 

Where a = F = ^m^^e 2 ' From tne Fq- (21) we can get t ne relation 
between the resonance amplitude and the shift frequency from the resonance 
point. 



Fig. 3 

The figure shows the relation between the square of resonance amplitude a and 
the shift frequency from the resonance point A, which is the solution of the Eq. 
(21). If the amplitude of noise F is small, the peak occurs at A = 0, as the 
amplitude of the noise increases, the shape of a changes gradually while it keeps 
the character that has the only one maximum, which truly shows the existence 
of the stochastic resonance between the myosin head and the thermal noise. 

After a sufficient time, the damping terms of the Eq.(16) go to zero, and only 
oscillational term ae 1 ^* - ^ survives. It is clear that these equations (Eq. (1)- 
Eq. (4)) do not have an energy source except the thermal noise term F(t). In 
the case, the energy of the myosin has the constant value which is determined 
by the amplitude a. And the average absorbed energy per unit time will be 
E(to) = 2fl 2 fima 2 . 



4 Numerical results 

We adopt a numerical method to solve the above equations. The parameters of 
the equations are shown in the Table. 1. 



Table. 1 
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The numerical results are shown in Fig. 4-Fig. 8. 



Fig. 4 

The figure shows the horizontal position of myosin head, the trace has irregular- 
ity and randomness like Brownian particles, but movement is a translacational 
one as a whole. 



Fig 5 

The solid line shows the vertical motion of the myosin head and the stretching 
motion of the rod is described with the dash line. From the figure we can know 
that the stretching vibration of the spring has irregularity and the collision be- 
tween the G-actin and the myosin head distorts their relative motion from a 
trigonometric function. 



Fig 6 

The dash line shows the motion of the filament and the solid line shows the 
position of the G-actin. From the figure we can know there is a relative sliding 
motion between the filament and the G-actin. The filament can move to one 
direction while the G-actin moves to opposite direction. 



Fig 7 

The figure gives the relative between the horizontal movement of myosin fila- 
ment and the initial angle. From the figure we can know that the motor can 
move fastest when the initial angle is nearly equal to 45° . But if the initial 
angle is too small or too big, the motor moves slowly, even cannot move. So the 
initial angle of the inclined rod is important for SIRM to move to one direction. 



Fig 8 

From the figure we can know that the shape of the G-actin is important for the 
system to move. If the G-actin if flat (R = 0) the SIRM can not move, on the 
other hand, if the radius of the G-actin is too big (R — 40). The system can 
not move, either. The motor can move at the highest speed at R — 20. 



5 Summary and conclusion 

Based on SIRM we present a sliding mechanism for the actin myosin system 
and set up the dynamic equations for the model. Because equations (Eq. (1) 
- Eq. (8) ) are so complex and highly nonlinear that we cannot obtain a strict 
analytical solution, we adopt a numerical method to solve these two dimensional 
Langevin equations. We discuss the motion of the myosin head and find there 
is a stochastic resonance between the head and the noise. Our model propose a 
microscopic mechanism for the actin myosin system: firstly, the thermal noise 
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interacts with the myosin head, and the resonance occurs between the myosin 
head and the noise, then it collides with a G-actin and obliquely kicks the G- 
actin sphere, because the direction of vibration is inclined against the line of the 
actin fibers, myosin molecules obtain the propellant force along the direction of 
the fiber and the myosin can move to one direction. 

From the numerical results we can know there is a relative sliding motion 
between the filament and the G-actin. The filament can move to one direc- 
tion while the G-actin moves to opposite direction. When the system absorbs 
the energy from the thermal noise constantly through stochastic resonance, the 
intermolecular potential and the inclined rod make the filament move to one 
direction. The system can convert the random noise to one directional motion. 
SIRM is thermally open to the outer surroundings and it has a outer heat source 
of ATP. The heat consumed by SIRM is directly supplied from the surround- 
ings and the energy provided to the surroundings comes from the outer source 
or reservior through a general energetic flow. From a macroscopic point of view, 
SIRM can move by getting heat from its surroundings. 

Noise in dynamical systems is usually considered a nuisance. But in certain 
nonlinear systems including electronic circuits and biological sensory apparatus 
the presence of noise can in fact enhance the detection of the weak signals. From 
the above results, we point out the possibility that the actin myosin system can 
use thermal noise around for the movement in similar way. This viewpoint is 
completely new as the traditional viewpoint thinks the water disturbs the mo- 
tion as viscosity. 
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6 Figure captions 

Table. 1. The parameters of the equations, a, (3, 77, 7 are viscous constants, 
M, m, M g are the mass of myosin bundle and myosin head, respectively, L is 
nature length of the spring, p, q, U are the parameters of the potential, Kg, Ki 
are tangent constant and radial constant of the spring. 

Fig. 1. Actin myosin system as an element of the muscle shown in the bottom 
of the figure. The structure of muscle is also shown from the upper level to 
lower level (see ref.[21]). 

Fig. 2. The model for actin myosin system. The coordinates of the head of 
myosin and the end of myosin connected to the filament are {x, y) and (x 2 , y 2 ), 
respectively. The body of the myosin is assumed to be represented as a spring. 
Fig. 3. The figure shows the solution of Eq. (21) and give the relation be- 

o 2 

tween the square of resonance amplitude a = ^ and the shift frequency A 
from resonance point. From the up to the bottom, different lines correspond to 
F = 32m4g e 2 : 20 ' 15 ' 10 ' 5 ' 1 ' respectively. 

Fig. 4. The figure shows the horizontal motion of the myosin head (see Eq. 
(I))- 

Fig. 5. The solid line shows the vertical movement of the myosin head 
and the dash line presents the stretching displacement s of the rod and s = 
^/(x-x 2 y + {y-y 2 y-L (see Eq. (2)). 

Fig. 6. The dash line gives the movement of the filament (x 2 ) and the solid 
line shows the position (x g ) of the G-actin (see Eq. (3) and Eq. (4)). 
Fig. 7. The curve shows the horizontal displacement of the myosin filament at 
different initial angle 6*o (t = QOunittime) . 

Fig. 8. The curves give the relation between the radius of G-actin and the 
displacement of filament in given time, different lines correspond to different 
radius R shown in the figure. 
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Table. 1. The Parameters of the Equations. 



1 unit time 


10- 1L s 


a,P,V,1 


1 


U 


20000 


1 unit length 


i(r 9 m 


p,q 


100 


m 


1000 


1 unit mass 


2.5 x 10" 2i kg 


L 


100 


M 


3000 


Kg 


500000 


Kt 


500 


M g 


10000 



11 



This figure "Figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0306192vl 



This figure "Fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0306192vl 
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